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Abstract 

We examine the Constrained Minimal Supersymmetric Standard Model (CMSSM) 
with an additional requirement of strict b-r unification in the region of small tan/3. 
We find that the parameter space becomes completely limited below about 1 TeV by 
physical constraints alone, without a fine-tuning constraint. We study the resulting 
phcnomenological consequences, and point out several ways of falsifying the adopted 
b-r unification assumption. We also comment on the effect of a constraint from the 
non-observation of proton decay. 
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1 Introduction 



A recent resurgence of strong interest in supersymmetry (SUSY) has led to a number of 
attempts at exploring the physical spectra and phenomenological consequences associated 
with the Minimal Supersymmetric Standard Model (MSSM) in the context of Grand Unified 
Theories (GUTs). This renewed interest is due primarily to the realization that gauge 
coupling unification within the Standard Model (SM) does not occur at any scale as one 
would expect from GUTs such as SU(5). On the other hand, within the MSSM such 
unification is found to be possible. Early studies H have been followed by a series of 
increasingly elaborate, and increasingly precise, analyses which have built complete SUSY 
spectra consistent with the unification assumption (see ]2| and references therein). These 
studies have mostly used the well-motivated supergravity (SUGRA) assumptions which 
suggest equating many of the unknown soft SUSY-breaking mass terms at the GUT scale. 
Within this framework, very complete studies can be done covering the entire range of 
possible SUSY masses, and specific, testable predictions can be made. 

In a previous work Q, we have examined the MSSM under a number of general assump- 
tions about the unification of the gauge couplings and masses, independent of the choice of 
gauge unification group.Q Within this context, a number of predictions, bounds, and signals 
were deduced and studied. A specific choice of GUT model and/or further assumptions or 
constraints could only serve to sharpen these predictions. 

In this letter we consider one further aspect of unification: the apparent unification of 
the bottom quark and tau lepton Yukawa couplings at the GUT scale, as would be expected 
in a GUT containing minimal SU(5) Yukawa interactions. Semi-analytic studies completed 
recently by several groups ||, ||], || have found, however, that the MSSM does not in 
general produce this b-T mass unification except in small and well-defined regions of the 
parameter space of mt and tan/3. Specifically, it has been realized that, up to GUT-scale 
threshold corrections, b-r mass unification can only occur if the top Yukawa coupling is at 
or near its infrared pseudo-fixed point. 

In Ref. [0], we examined the size of the GUT-scale corrections necessary in order to 
remove the strong constraints on tan/3 and mt- We found that 0(10%) corrections both in 
the gauge and Yukawa unification were more than adequate for allowing b-T mass unification 
over very wide ranges of values for tan/3 and mt- Nonetheless, examinations of "typical" 
GUT threshold corrections || have yielded corrections too small to significantly alter the 
relation between tan/3 and mt, so a detailed exploration of the SUSY parameter space 
indicated by b-T mass unification seems well-motivated. Our goal is to derive testable 
experimental consequences of this assumption, and to point out how b-r mass unification 
within the MSSM can be falsified in a number of ways. Finally, we will comment on the 
effects of a constraint from the non-observation of the proton decay. 

2 The Pseudo-Fixed Point Solutions 

It was recognized several years ago that one of the strengths of the MSSM over the SM 
was that in the MSSM the b- and r- Yukawa running couplings meet at roughly the same 

1 We only require that the choice of unification group lead to sin 2 9 VJ (Mx) = § which also holds in many 
phenomenologically viable superstring-derived models. 
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mass scale at which the unification of the gauge couplings takes place ||. As both the 
experimental data and the theoretical calculations became more precise, it became clear that 
gauge coupling unification within the MSSM occurs over the entire range of theoretically 
acceptable SUSY mass scales. At the same time, however, Yukawa unification within the 
MSSM is not so trivial. For most choices of input parameters (e.g., tan/3 and mt), the r- 
Yukawa coupling is as much as 0(20%) larger than that of the 6-quark at the gauge coupling 
unification scale. And because the slopes of the bottom and tau Yukawas are typically flat 
at large scales, the Yukawa couplings can "unify" (i.e., cross) at a mass scale many orders 
of magnitude smaller than the GUT scale, even though their GUT-scale values may differ 
by only 10-20%. 

The exception to this general trend, however, occurs when either (i) the top Yukawa 
coupling, or (ii) the b- and r- Yukawa couplings, are at or very near their infrared pseudo- 
fixed point values. That pseudo-fixed point value is the value of the Yukawa coupling (at 
the electroweak scale) which produces a Landau pole precisely at the GUT scale; that is, 
values of the coupling greater than the fixed point value will become non-perturbative at 
scales below Mx when run up using its renormalization group equations (RGEs). 

The values for the top, bottom, and tau Yukawa couplings corresponding to the pseudo- 
fixed point can be precisely determined. For the case at hand, one finds @, f§ that there 
are two conditions for sitting on or near one of the pseudo-fixed points (that is, for finding 
b-T Yukawa unification) either one of which must be satisfied. Either 

M t ~ (200 GeV) sin 3 (1) 

or 

50<tan/3<70. (2) 

The first condition corresponds to the top Yukawa pseudo-fixed point, while the second is 
the bottom-tau Yukawa pseudo-fixed point. Also, because the difference between the top 
quark pole mass and running (MS or DR) mass can be as much as 10 GeV in the region of 
interest, one must be careful to specify to which top quark mass one is referring. We will use 
TYiq to specify a quark running mass (or a general definition when the choice is irrelevant) 
and M q its pole mass. 

In Fig. [l] we have shown the region in the (Mt, tan (5) plane consistent with b-T mass 
unification. The width of the region is due to the 3a uncertainty in the 6-quark mass (using 
M b = 4.9 ± 0.07 GeV [§) and the la uncertainty in a s (m z ) (a s (m z ) = 0.120 ± 0.007 @); 
in this figure none of the width comes from GUT threshold corrections to the Yukawa 
unification. 

From Fig. [l] it is evident that there are three primary regions of interest in this plane. 
The first is the region (labelled I) of large tan/3, over all Mt, where the b- and r- Yukawa 
couplings reach their pseudo-fixed points. Some part of this region, corresponding also 
to large Mt, is of particular interest to those studying S'O(IO) unification with a minimal 
Yukawa sector. There one can obtain the GUT-scale prediction m& = m T = mt, which 
when renormalized at the electroweak scale yields tan 8 ~ — ~ 50 ~ 70. Because of that 
relation, this region deserves consideration and some studies have been carried out [||, TO, Q|. 
However, certain difficulties invariably arise in considering the MSSM with very large tan/3. 
Three in particular stand out. 
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Figure 1: The (M t ,tan/3) plane showing the region consistent with b—r unification. The width of 
the region is due to the 3<r uncertainty in Mb and the la uncertainty in a s (see text). The three 
regions of unification are all visible: (I) tan/3 ~ 50 ~ 70, (II) 190 GeV < M t < 210 GeV, and (III) 
tan/3 ~ 1. 



First, it has been argued that within the MSSM (with two Higgs doublets) large tan (3 
is unnatural [12]. Specifically, one finds that large hierarchies, only some of which can 
be protected by additional imposed symmetries, arise among the mass parameters of the 
Higgs potential. Second, one has particular difficulty with the 1-loop corrections to the 
6-quark mass, which are proportional to tan/3 and can easily be of O(50%) unless new 
symmetries are imposed in order to control them KM. This issue has been more recently 



investigated in Ref. [11] in the context of an 5O(10) GUT model and it has been argued that 
requiring dynamical electroweak symmetry breaking makes those corrections well defined 
though generally unsuppressed. Finally, models with large tan (3 tend to produce very small 
branching ratios for the FCNC process b — > S7 due to their suppressed charged Higgs 
contributions, perhaps inconsistent with recent CLEO data. Therefore, we choose to put 
off any further consideration of this region for now. 

The second region (labelled II in Fig. [I]) leading to b-r mass unification is found at large 
M t between about 190 and 210 GeV for almost all values of tan (3 (between about 2 and 60). 
This region of parameter space is strongly disfavored by the 2a LEP limit mt < 180 GeV 
for m h o < 120 GeV ]^q| . Further, if FNAL does see top quark events, one could exclude 
top quark masses in this region due to the small cross section expected for such large mt- 
Therefore, we will not consider this region in the present study. 

The third and final region (labelled III in Fig. |) is at low tan/3, with M t < 190 GeV, 
where there is almost a one-to-one correspondence between tan/3 and Mf. This range is 
perhaps of more immediate interest since it corresponds to the top mass range to be covered 
by the Tevatron. And, having disfavored the previous two regions, we are led to consider in 
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detail this region of tan (3 near one, for 155 GeV < Mt < 185 GeV. As we will show below, 
this scenario leads to very specific and falsifiable predictions which could rule it out. 

One should remark on the size of the radiative corrections to the 6-quark mass in the low 
tan (3 regime. We find for all solutions in this third region that the (leading) gluino-induced 
corrections to the bottom quark mass are always small, #m{,/m& < 2%; the higgsino-induced 
corrections are approximately four times smaller still. These corrections are far too small 
to alter our results. We also find that the sizes of these corrections depend only very weakly 
on the SUSY mass parameters mo, an d in fully consistent solutions (see following 

section), and so remain small even at large SUSY mass scales. 

3 The Constrained Minimal Supersymmetric 
Standard Model (CMSSM) 

In Ref. 0], we described the construction of what we termed the Constrained MSSM (or 
CMSSM), built by relating the MSSM soft-breaking terms through the minimal SUGRA 
assumptions and then constraining these solutions as summarized below. Here we briefly 
summarize the basic points of that construction. 

For each choice of M t , mo, mi/ 2 , Ab and sgn//o (tan/3 is now determined through the 
requirement of b-r unification as described earlier) we find a solution in the CMSSM; each 
solution is a complete set of values for a s and the masses of the Higgs bosons and all the 
superpartners. (The exact procedure for building such complete solutions is summarized 
in Ref. j|].) If the Higgs potential at the electroweak scale does not admit electroweak 
symmetry breaking (EWSB) consistent with the SM, that choice of input parameters is 
discarded. Further constraints are then applied. 

Besides requiring that EWSB occur, we demand that all physical mass-squares remain 
positive. We demand that all solutions be unobservable by current direct experimental 
searches, including the requirement that the solutions provide a BR(6 — ► S7) consistent 
with CLEO data. Furthermore, we calculate the relic density of the lightest SUSY particle 
(LSP), demanding only that the LSP be neutral, and, from limits on the age of the universe 
of 10 billion years, we demand that QLSpho < 1. Those solutions which finally remain after 
all these cuts comprise the allowed parameter space of the CMSSM. 

4 Results 

We have examined solutions for three representative top quark (pole) masses of Mt =155, 
170, and 185 GeV in order to study the range that seems to be indicated by the LEP analysis 
and early CDF results. In building the parameter space of solutions, we have varied mo 
and mi/2 from 20 GeV to the TeV scale logarithmically; no upper bound is set by hand on 
mo or mi/2) their upper bounds coming from the constraints on the solutions which define 
the CMSSM. The lower bound of mo was not taken to zero, but values in the lower region 
(mo — O(20 GeV)) are for all purposes phenomenologically identical to each other as well 
as to No-Scale models, since experimental bounds force ^ 80 GeV. The value of Aq 
was restricted to the range —2.5 < Aq/tuq < 2.5 in order to avoid possible color-breaking 
minima of the full scalar potential. The final free parameter, sgn/Lio, was allowed to take 
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both possible values. 

4.1 Allowed Parameter Space and Mass Spectra 

One of the remarkable consequences of considering this low tan (3 limit with a complete 
analysis is the existence of upper bounds on the mass parameters of the MSSM without 
resorting to an imposed arbitrary fine-tuning condition. In particular the combination of 
various mass bounds from direct searches, the age of the universe constraint, and the re- 
quirement that the LSP be electrically neutral combine to put strict upper (and lower) 
bounds on mo and mi/2, and therefore on the mass spectrum of the MSSM. These bounds 
are entirely due to the physical constraints of the CMSSM. As we pointed out previously Q, 
for small tan/3 < 2 and/or large mj > 170 GeV, the parameters m\/2 and m,Q (and there- 
fore the whole SUSY spectrum) are completely constrained from above in the O(lTeV) 
range. The case considered here falls into that category. These absolute upper bounds are 
m^-dependent and are usually somewhat weaker than those imposed by simply requiring all 
SUSY masses below about 1 TeV, or than those which the fine-tuning constraint we used in 
Ref. [|j would have permitted. For this reason, in examining some phenomenological appli- 
cations of the solutions in the CMSSM, we will place an additional fine-tuning constraint. 
We do so in order to ensure that phenomenological results of this study are "realistic" ; that 
is, although consistent solutions may exist with large SUSY mass scales, we wish to exclude 
these from phenomenological consideration on the basis that they reintroduce too much 
fine-tuning into the physics. (Our precise definition of fine-tuning is discussed in Ref. Q. 
We note that the definition that we use diverges at tan/3 = 1; the values of tan/3 that we 
consider here are far enough from unity so that this effect is not significant, and becomes 
irrelevant as tan/3 increases with increasing mj.) 

In Fig. |2] we have shown the allowed range of parameter space in the (m 1 ^, mo) plane for 
M t = 155 (Fig. |2|a), 170 (Fig. ^b), and 185 GeV (Fig. ^c) without fine-tuning constraints. 
All three graphs exhibit many similarities which are general features of the CMSSM for 
all tan/3 Q. In both cases the region of large myz 3> mo is excluded by demanding a 
neutral LSP. This uniquely selects the lightest neutralino % = Xi (mostly a bino-like state) 
for which we calculate the relic abundance. Then we find that large m-o are excluded by 
requiring ^ x /iq < 1- 

For the Mt = 155 GeV case, we find that Yukawa unification allows tan (3 to take on 
values only in the range 1.1 < tan/3 < 1.4. This, coupled with the requirements of the 
CMSSM, restricts 80 GeV < m 1/2 < 940 GeV and m < 350 GeV. When we apply the fine- 
tuning constraint / < 50, we find the approximate bound shown as a solid line in Fig. ^a. For 
this subset of solutions, we find new upper bounds on the parameters of the model. In par- 
ticular, the fine-tuning constraint yields tan/3 < 1.2, myz ^5 180 GeV and mo < 210 GeV. 
Notice also the strong lower bounds placed on the parameter space. This bound comes 
from two sources: for [i > 0, solutions with small m^2,mo are ruled out by the Higgs mass 
bound; for fj, < 0, it is the t\ mass bound that rules out the same approximate region. This 
effect is strongly diminished with increasing tan/3 (and therefore Mt). 

The Mt = 185 GeV case in Fig. [2|c exhibits one additional interesting feature. For a 
relatively narrow range of my 2 — 100 GeV solutions with much larger mo are allowed. 
The source of this behavior is the Z-pole enhanced neutralino pair annihilation into a pair 
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Figure 2: The regions of the (mi/2, mo) plane consistent with low tan/3 b-r mass unification, given 
all the constraints of the CMSSM. For (a) M t = 155 GeV, (b) M t = 170 GeV and (c) M t = 185 GeV. 
Each dot represents one solution belonging to the CMSSM. Because mo and mi/ 2 are discrete inputs 
in this approach, the points have been "smeared" to show variations in density. The regions with 
no solutions are labelled as to which constraints forbid solutions there: age of universe bound (A), 
neutral LSP requirement (L), non-tachyonic stop (T); Higgs mass bound (H), and chargino mass 
bound (C). The solid curves in each figure are the approximate cutoffs dictated by our choice of 
fine-tuning constraint. 



of ordinary fermions. (For m^ 2 — 100 GeV, m x ~ 0.45 m^2 — mz/2.) This exchange 
vanishes in the limit tan/3 = 1 and is therefore negligible in Fig. ^a where tan/3 < 1.4. 
For Mt = 185 GeV, we find 1.8 % tan/3 < 3.1 and the effect of the Z-exchange becomes 
important. Our analysis provides an overall envelope of 100 GeV < my 2 % 1.1 TeV and 
niQ < 600 GeV. Once again the fine-tuning constraint tends to lower the upper bounds 
of the various model parameters. In particular, we find tan/3 < 2.4, m 1 / 2 ^5 290 GeV and 
mo 420 GeV after applying the fine-tuning constraint. The corresponding fine-tuning 
bound is shown solid line in Fig. 0c. 

The intermediate case for Mt = 170 GeV in Fig. does demonstrate some hint of 
the Z-pole effect found at the larger tan/3 associated with Mt = 185 GeV. Here we 
find 1.4 < tan/3 < 1.9, 90 GeV < m 1/2 < 940 GeV, and m < 500 GeV without a fine-tuning 
constraint. Upper bounds with the fine-tuning constraint are modified to be tan/3 < 1.5, 
m 1/2 < 260 GeV and m < 300 GeV. 

It is worth noting that the region itlq — fnz (and small tan/3) is favored by 

the non-observation of proton decay in S'C/(5)-based GUTs with minimal Higgs sector [14|. 
In this case the neutralino relic abundance is reduced by the Z- and /i°-pole effects. In 
this region the LSP still remains mostly bino-like, although with somewhat smaller bino 
component ( > 90%). However, predictions for proton decay can be suppressed with more 
complicated Higgs sectors (see, e.g., Ref. Ilql ) and, since we do not assume any specific 
GUT model here, we will also not use this constraint to limit the parameter space of the 
CMSSM. We also note that our numerical routine for the relic abundance is not designed 
to properly calculate ^x^o ^ n ^ e vicinity of a pole (that is, within about 10 GeV of the 
pole) and therefore the values of ^ x /iq in this regions are only indicative. 

One should note for all bounds throughout this study that the exact values depend on 
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Table 1: The bounds of the masses and parameters of the MSSM under the constraint of low tan/3 
Yukawa unification for M t — 155, 170, and 185 GeV. For each top quark mass the lower bound 
(labelled lower), the absolute upper bound with no fine-tuning constraint [upper) and the upper 
bound with the fine-tuning constraint (FT) are shown. All masses are in GeV. The general masses 
q and I l represent the bounds on all squarks and LH-sleptons excluding the third generation. The 
bounds on all first and second generation RH-sleptons are essentially those of the ti for low tan/3. 
Note that there is sensitivity to the grid of values for (too, m-i/2, A)) that we have chosen, and 
therefore some uncertainty in the exact bounds. 

our numerical sampling of the original input parameter space and so should be considered 
with appropriate errors. In particular, upper (lower) bounds on mx/2 could be modified 
upwards (downwards) by as much as 12% with a smaller sampling grid; bounds on mo 
could likewise be increased (decreased) by as much as 20%. 

Because such strict bounds exist for these cases (with and without fine-tuning), we can 
place bounds on the CMSSM mass spectra. In Table [l] the numerical bounds on a variety of 
important quantities are shown for all three top masses, with and without the fine-tuning 
constraint. 

It is significant that much of the region of low m^ for the Mt = 155 GeV case has been 
excluded on the basis of the non-discovery of the Higgs boson in direct searches. For all 
solutions in the superunified MSSM, the lightest Higgs scalar (h°) has essentially identical 
properties to those of its SM counterpart; that is sin 2 (/3 — a) ~ 1 always. Therefore 
mass bounds on the SM Higgs apply equally well to the h°. (We conservatively require 
m h o > 60 GeV.) However, h° receives large one- loop radiative corrections which can increase 
its mass by 0(mz)- Therefore, limits on the solution space based on Higgs non-observation 
must include the full one- loop radiative corrections, which in turn requires a full calculation 
of the SUSY mass spectra, particularly the mass eigenstates of the top squarks. 

Because the top mass dependence of these radiative corrections is so strong, one can 
clearly understand why the lower bound on mi/2 due to the Higgs mass decreases with 
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Figure 3: The number of solutions with a given m^o versus the value of m^o, for M t = 170 GeV 
(dashed) and 185 GeV (solid). Only the solutions with small fine-tuning are plotted. 

increasing mt ; for larger mt, one gets larger m^o from radiative corrections, so the exper- 
imental bound on m^o rules out less of parameter space. In Fig. || we demonstrate this 
behavior by plotting the number of solutions leading to a given m h o versus m h o, both for 
Mt = 170 GeV (dashed) and 185 GeV (solid), subject to our fine-tuning constraint. 

There are also upper bounds on m^o, and several groups have recently examined these 
bounds within this approach (l^, [R], 17|. However, because in the tan/3 — > 1 limit m 2 h0 — > 



at tree-level, these bounds are highly dependent on the size of the radiative corrections to 
m^ . And because these corrections increase quartically with mt and logarithmically with 
the SUSY masses, their size is highly dependent on one's assumptions about how to cut off 
the allowed MSSM parameter space, that is, how one defines what "too much fine-tuning" 
means. 

For M t = 155, 170, and 185 GeV, we find m h o to be less than 81, 99, and 118 GeV 
respectively, with the fine-tuning condition / < 50 imposed. Without this condition, the 
upper bounds increase, as the effective SUSY mass scale is somewhat increased (compare 
Table |l|). Though we believe that it is the first set of bounds that should be taken as more 
indicative of our expectations since they more fully contain theoretical prejudices which 
apply to SUSY, one should bear in mind the strong dependence of these bounds on the 
choice of fine-tuning condition. This caveat, however, does not apply to the maximum 
values as given in Table [l], where no fine-tuning condition at all was used to bound the 
parameter space. One should also note that calculations of the two-loop corrections to m^o 



show a net decrease of m h o below its one- loop value [18], and therefore will not disrupt our 
bounds. 

Detection of the other Higgs bosons (H°, A , and H^) cannot be accomplished at 
LEP II nor at the proposed NLC500 0. We find for 155 GeV < M t < 185 GeV that 
m H o ~ m^o ~ m H ± > 260 GeV, outside the range of either machine. We believe that the 
detection of the heavier Higgs scalar H° might be possible at the LHC somewhat beyond 
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the asserted region m H o < 2Mt pQ ]. However, all solutions (with small fine-tuning) do have 
at least one SUSY particle that is detectable at the LHC in addition to the h°, even if the 
heavy Higgs bosons are not. 

Can LEP II find any sparticles? We find that under the assumption of low tan (3 b-r 
unification, detection of states other than the lightest Higgs may be possible, but for a few 
solutions only. Searches for light SUSY particles should concentrate on the lighter stop (tx), 
the lighter stau the lighter chargino (xt) and the second lightest neutralino (x^j where 
m o ~ m ±). Within the range of Mt considered, we find solutions with masses for these 
particles all the way down to their current experimental limits. In particular, for such light 
stops there is a large mixing of the right-handed and left-handed interaction eigenstates, so 
the simple approximation that t\ ~ £_r does not hold. Further studies of the detectability 
of the MSSM under similar constraints and assumptions have been done in Ref. |i"7| , fzTJ . 

We finally note that the resulting ranges of m~ and scalar quark masses are largely 
above the reach of the Tevatron. Thus finding any of those particles well below the ranges 
indicated in Table || would rule out the b-r unification if tan j3 is close to one. 



4.2 BR(£> -> si) 

In the general superunified MSSM where Yukawa unification has not been required, the 
recent CLEO bounds on BR(6 — > 57) have the ability to rule out certain regions of parameter 
space and indicate a future ability to further constrain or discover SUSY through more 
precise measurements of BR(6 — >• S7) Q. However, one finds in the region of low tan/3 
Yukawa unification that almost all solutions consistent with all other requirements of the 
CMSSM naturally fall within the bounds of the CLEO data, and in particular, no solutions 



provide larger branching ratios than are allowed by the data. (We follow Ref. [22] for our 
calculations of the branching ratio.) But as the following analysis emphasizes, this is not in 
general due simply to the decoupling of the SUSY contributions; in fact, for solutions with 
low fine-tuning, the SUSY and SM contributions are often comparable in size. 

Nonetheless, what is particularly noteworthy in the low tan (3 and small fine-tuning 
limit is that the branching ratio is highly dependent on the sign of fi. In Fig. || we have 
histogrammed our calculated BR(6 — * 57) for Mf = 170 and 185 GeV solutions with small 
fine-tuning (we comment on relaxing the fine-tuning condition below). The central peaked 
region falls approximately at the SM prediction of the branching ratio, with larger Mt 
moving the peak (and the SM prediction) to larger values. Note however that the histogram 
yields two separate, non-overlapping regions. In both cases, the region of higher branching 
ratio is occupied only with solutions of \i > 0. Likewise the region of lower branching ratio 
is occupied only with solutions of jx < 0. Presuming that both theoretical and experimental 
uncertainties to b — > 57 ever become small enough, one may then have a method by which 
to differentiate the sign of \i through this process. 

Why the two separate regions? Of the non-SM contributions to the b — > 57 amplitude, 
the dominant contribution in this region of parameter space tends to come from the X\~t\ 
loop, with sign opposite to that of the SM W -t contribution. For all acceptable solutions 
of the CMSSM, one gets a lighter chargino that is almost pure and so only has couplings 
to ti,. Without mixing of the t's it is the tn that invariably comes out to be the lightest 
squark, which leads to a very small xt~ti contribution 
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Figure 4: Histogram of BR(b -> s-y) for (a) M t = 170 GeV and (b) M t = 185 GeV, for solutions 
corresponding to low fine-tuning. 



Now it is well known that the ti and tn can have large mixings proportional to m 2 (At + 
n/ tan (3) 2 . For small mixings the lighter stop (t\) is almost pure £r; as the mixing increases 
t\ gains a larger ti component and so the Xi~ti contribution becomes sizeable. For solutions 
in which tuq is not too much larger than mj/2, one finds that A% is driven negative at the 
electroweak scale regardless of its magnitude and sign at the unification scale. Therefore, if 
fi is also negative, large ti-tn mixings result and t\ has a sizeable ti component; if fi > 0, 
the A t and n contributions to the mixing tend to cancel, forcing the ti component of t\ to 
be very small. 

In the [i > case, then, the xt~ti coupling will be very small, allowing the SM con- 
tribution to easily dominate. In the \i < case, however, the coupling will be sizeable, 
cancelling much or all of the SM contribution. Therefore, the \x < case will give much 
smaller branching ratios than one would get from the \i > case, with the two regions 
separated near the SM prediction. 

What happens as we allow larger tan/3 or larger fine-tuning? Since the i^-is mixing 
goes as (At + fi/ tan/3) 2 , it is clear that as tan/? increases, At will come to dominate the 
mixing and the results will be relatively independent of sgn/xo- Similarly, the dependence 
on sgn disappears for solutions corresponding to larger fine-tuning because in this case (i 
typically becomes large and its contribution dominates over that of A t , producing mixings 
proportional only to (mtfi) 2 . However, as expected, for solutions corresponding to large 
fine-tuning, and therefore large masses in the loops, the supersymmetric contributions go 
to zero, leaving only the SM contribution and thus no sgn fio dependence anyway. 

4.3 The ^-parameter 

Of particular interest is the Higgs/higgsino mass parameter ji which does not break SUSY 
and therefore could in principle take values much larger than the soft SUSY-breaking param- 
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Figure 5: Scatter plot of M2 versus \i(m.z) for (a) M t = 155 GeV and (b) M t = 185 GeV for all 
solutions in the CMSSM without a fine-tuning constraint. Notice the "mo ^> nii/2" points present 
in (b) allowed because those solutions have LSP's falling in the pole of the Z-channel annihilations. 
The tree level calculation of the \1-M2 correlation is also shown (dashed line). 



eters. In this approach, however, its size is determined through the condition of electroweak 
symmetry breaking and comes out to be of the same order of magnitude as m.1/2 and mo, 
as has been discussed in detail in Ref. Q and many other places. Further potentially strong 
correlations can be derived by imposing additional constraints or assumptions, like the b-r 
unification discussed here. 

Working in the top Yukawa pseudo- fixed point limit, the authors of Ref. [23] found 
semi-analytic expressions indicating a strong correlation between fj, and m^j in the region 
of Yukawa unification. Though their results were only at tree level and so did not contain 
contributions from the one-loop effective potential, they are easily extended to include the 
leading correction from the t sector. We find that for m.1/2 > fnoj m z'- 



2 2 0.5 + 3.5 tan 2 /3 



15«2 

2vr 



mi 



m 



W , 



I 'l l % 



(3) 



The first term on the right in Eq. (g) is the tree level contribution only [23], while the 
second term represents the one-loop leading correction. 

We illustrate the behavior of fi as a function of M2 (M2 ~ 0.8mi/2) in Fig. [5] for our 
solutions with (a) Mt = 155 GeV and (b) Mt = 185 GeV. One can compare this to Fig. 31 
of Ref. Q where the (/j,M2) plane for general solutions in the CMSSM was displayed. 
There one does not see evidence for the strong correlation between ^ and M2 (or jti^) 
that one finds in the pseudo-fixed point limit. We have also plotted in Fig. [| a dashed line 
corresponding to the tree level calculation of fi as given by the first term in Eq. (||). 

From Fig. |B| one sees that the tree level expression for the n-mi^ correlation describes 
our solutions well until M2 ^ 400 GeV, where the slope rises due to the one-loop corrections. 
The one- loop effects are large enough so that the tree level calculation for \x at m^ 2 — 1 TeV 
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is 0(50%) larger than the actual one-loop value. 

One should note, however, that there is an ambiguity in the choice of scale at which 
one renormalizes the SUSY masses versus the scale at which one minimizes the potential, 
leading to uncertainties in the one-loop contributions that can be large. The problem stems 
from the fact that we minimize the potential at Q = mz, while the SUSY masses are 
renormalized at their thresholds. This leads to corrections of the one-loop potential that, 
though of two-loop order, can become significant [|l^, |24j]. Nonetheless, because the solutions 
with large one-loop contributions are disfavored by fine-tuning arguments, one can safely 
ignore this question. 

The only significant deviation from Eq. (J3J) occurs in Fig. ||b for some points at very 
low m-i/2- These points correspond to the mo S> m^ 2 points in Fig. ^ where the relic 
density was suppressed due to the presence of the Z- and /i°-poles. These points can be 
missed in semi-analytic approaches without calculating the neutralino relic density. It is 



also this region which S\JSYSU(5) proton lifetime calculations like those of Ref. [14] select. 
For these points we are far from the limit in which Eq. (|3|) holds and one finds instead for 
rriQ ^> mx/2,Tnz another tree- level correlation between fi and mo |23| which fits our solutions 
well (the effects of the one- loop effective potential here are very small): 

o 9 1 + 0.5 tan 2 3 

tan^-1 • (4) 

Finally, one also sees from Fig. |5| that one can put absolute upper and lower bounds on 
within this framework, without regard to M2. We have included these Mt-dependent 
bounds in Table [l] and Fig. |5[ 



4.4 Neutralino Relic Abundance and Dark Matter 

We have already emphasized the crucial role played by cosmological constraints in deriving 
upper bounds of 0(1 TeV) on the supersymmetric mass parameters, without having to 
impose a somewhat arbitrary constraint of fine-tuning. The main ingredients that lead to 
such upper bounds in the CMSSM are: (i) the lightest neutralino \ comes out to be the 
only possible LSP which is neutral; and (ii) it turns out to be predominantly bino-like. 

Even after rejecting solutions with charged LSPs as dark matter (DM) candidates (in 
our case it is t\ in the region m^ 2 ^> mo), the sneutrino could still have come out to be the 
(neutral) LSP. However, just as in the more general case without Yukawa coupling unifica- 
tion pi, we never find that to be the case after we apply experimental limits. Furthermore, 
the neutralino LSP comes out to be almost pure bino (p~ > 95% except for the region of 
mo S> mi/2 where > 90%) because of the radiative EWSB requirement which effectively 
leads to > M2 (see Fig. ||). 

Finally, it should be appreciated (even if it has been known for some time) that a 
neutralino relic abundance ^x^o c l° se to unity corresponds to sfermion masses in the range 
of a few hundred GeV, which is a natural mass scale in the MSSM with softly broken SUSY 
and radiative EWSB. This fact makes the neutralino an excellent candidate for the dark 
matter in the universe. 

Since in this approach all the masses and couplings are determined in terms of just a 
few basic parameters, we can also reliably calculate ^l x ^o as a function of those same input 
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Figure 6: The regions of the (tox/2, tuq) plane consistent with low tan (3 b-r mass unification, given all 
the constraints of the CMSSM, for Mt = 170 GeV, Aq/ttiq = and /i < 0. Solutions outside the thick 
solid lines are excluded: on the left (small mi/ 2) by the chargino mass bound (C) m x ± > 47 GeV 
and by tachyonic Vs (T); on the right (large m 1 / 2 3> mo) by charged LSP (L); and from above by 
the age of the universe, i.e. fi x /io < 1 (A). We also indicate the sub-regions selected by either the 
hypothesis of cold dark matter (0.25 < ^x^o ~ 0-5, between thin solid lines) or the one of mixed 
dark matter (0.16 < ^x^o ~ 0.33, between thin dashed lines). 



parameters. We include all the final states in calculating the neutralino pair annihilation 
even though the dominant contribution in most of the parameter space (away from the 
poles) comes from the exchange of the (lightest) sfermions. 

We reject those solutions for which Q x h,Q > 1 as corresponding to the universe being too 
young (less than about 10 billion years). This requirement alone appears to be extremely 
powerful, excluding values of m\/2 and tuq bigger than roughly 1 TeV and thus making 
low-energy supersymmetry a unique outcome of the simplest SUSY grand-unification as- 
sumptions. This is illustrated in Figs. ^ and |(| 

Furthermore, there is growing evidence for the existence of dark matter in the universe. 
While its amount and nature remain unclear, one of the most favored scenarios has been 
a flat universe (£1 = 1) with most of its matter (about 95%) contributed by DM. In one 
popular scenario the neutralino would be the predominant component of such (cold) DM 
in which case its relic abundance would be expected to be 

0.25 < Q x hl < 0.5. (CDM) (5) 

More recently (after COBE), a mixed CDM+HDM picture (MDM) has gained more atten- 
tion as apparently fitting the astrophysical data better than the pure CDM model. In the 
mixed scenario one assumes about 30% HDM (like light neutrinos with m v ~ 6eV) and 
about 65% CDM (bino-like neutralino), with baryons contributing the remaining 5% of the 
DM. In this case the favored range for is approximately given by 

0.16 < n x hl < 0.33. (MDM) (6) 
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Mass Limits 
(GeV) 


MDM 


CDM 


MDM/CDM 
+ FT 


lower 


upper 


lower 


upper 


m-i/2 


90 


520 


90 


660 


230 


m 


55 


245 


85 


245 


245 


\^( m z)\ 


270 


940 


280 


1120 


560 


M 2 


80 


430 


80 


545 


190 


h° 


62 


113 


62 


118 


98 




370 


1270 


380 


1530 


770 


*x 


95 


830 


82 


1020 


430 


n 


94 


OCA 

250 


110 


01 n 

31U 


250 


ii 


105 


380 


120 


480 


260 


q 


255 


1110 


265 


1350 


550 


9 


250 


1200 


250 


1480 


570 


X-LSP 


25 


225 


25 


290 


100 




48 


435 


48 


550 


200 


tan/3 


1.4 


1.7 


1.4 


1.8 


1.5 


a s (mz) 


0.123 


0.128 


0.123 


0.128 


0.125^ 



Lower bound with FT constraint. 



Table 2: The bounds on the mass parameters, tan j3, and a s of the MSSM under the extra constraints 
imposed by cold dark matter (CDM) and mixed dark matter (MDM) scenarios, for Mt = 170. The 
last column (MDM/CDM + FT) gives the upper bound when either the CDM or MDM scenario is 
assumed and the requirement of low fine-tuning (/ < 50) is additionally imposed. Uncertainties in 
the values are discussed in text and in caption of |l|. 



While neither of these DM scenarios is free from problems, it is nevertheless interesting 
to point out which regions of the parameter space of the CMSSM they select. This is 
illustrated in Fig. ^ and in Table |2[ We see in Fig. H that requiring either (|5|) or (^) results 
in selecting only relatively narrow bands in the (mi/2,Tno) plane. Their shape and location 
vary with other parameters but typically correspond to both and ttiq in the range of a 
few hundred GeV, independent of the choice of Aq and sgn^o- The resulting mass ranges 
are presented in Table ^ for Mt = 170 GeV. They should be compared with those listed 
in Table [l] for the same Mt to appreciate how much more limited the mass ranges become 
after the MDM/CDM assumption is made. It is clear that, with the exception of the light 
Higgs h°, the mass spectra consistent with either CDM or MDM are typically beyond the 
current experimental reach. Conversely, a discovery of a slepton at LEP II, or a squark 
(other than the stop) or gluino at the Tevatron well below the limits given in Table g, while 
providing unquestionable evidence for super symmetry, would at the same time indicate clear 
deficiency in the neutralino relic abundance [25| below the expected ranges of @ or (||), in 
the scenario with b-r unification and small tan j3. 

Finally, as we have noted previously, the tiny region mo S> Tn\/2 — mz, corresponding 
to the Z- and /i°-pole annihilation of the LSPs, is the region favored by the non-observation 
of proton decay in SU (5)-type models with minimal Higgs sectors [ff3. 
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5 Conclusions 



The predictability of the CMSSM becomes significantly enhanced by an additional assump- 
tion of b-r unification at least in the region of small tan j3 which we have studied here. Our 
main conclusions can be summarized as follows: 

1. The parameter space of the CMSSM is now completely limited both from below 
(by experimental constraints) and from above (mo % 500 GeV and m^j % 1 TeV) by 
robust cosmological constraints, without having to invoke the a fine-tuning constraint. 
This is a specific case of a more general property of the CMSSM: for either small tan (3 
(close to one) or for large M t > 170 GeV the parameter space is always (i.e., for any 
choice of other parameters) constrained from above broadly within the 1 TeV mass 
range. Both of these cases are selected by the requirement of low tan (3 b-T unification. 

2. The resulting sparticle mass spectra are highly constrained and correlated. All the 
colored sparticles (except for the lighter stop) are typically very heavy, and so are 
the heavier Higgs bosons (^4°, H°, H^). The lightest neutralino is the LSP, and it 
is invariably predominantly bino-like. Also, m x ± — m^o ~ 2m x . The resulting mass 
ranges, with and without imposing the additional fine-tuning constraint, have been 
listed in Table [I]. It is clear that one can make a number of predictions which can 
falsify the specific scenario considered here. It would be ruled out, for example, if Mt 
came out to be about 170 GeV and the gluino or a squark (other than a stop) were 
discovered at Fermilab well below 200 GeV; similarly, the sleptons cannot be much 
lighter than 65 GeV. 

3. Additional stringent constraints are provided by requiring the neutralino LSP to pro- 
vide most of presumed dark matter in the flat universe, as we can see from Table |2|. 
Again, we find that the lower bound on the slepton masses (including the stau) are 
beyond the reach of LEP II, while the squarks and the gluino could possibly be dis- 
covered with the upgraded Tevatron in a limited region of parameter space. Thus, 
finding such sparticles with masses well below the ranges given in Table | would pro- 
vide us with important information about the status of the neutralino as the dominant 
component of DM in the universe. 

4. The predictions for BR(6 — > 57) in this scenario fall almost completely within the 
range favored recently by CLEO, and near to the SM prediction. Furthermore, for 
solutions with light spectra a sharp dependence arises in the prediction of BR(6 —* 57) 
on sgn/uo- However, both theoretical and experimental uncertainties must be reduced 
before one will be able to constrain a large portion of the parameter space or determine 
sgn /io through this signal. 

5. In this restrictive scenario with imposed radiative electroweak symmetry breaking an 
additional correlation between [i and my 2 arises (see Eqs. (^)-(Q)) which may be 
helpful in a limited way in various phenomenological studies. 

Finally, which of the properties of the CMSSM sparticle spectra and predictions are spe- 
cific to the b-r unification assumption? Essentially, the crucial ingredient is the requirement 
that tan (5 be close to one. In this case the tree-level contribution to m h o is negligible and 
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h° is light enough to exclude large regions of the parameter space for smaller M t . Also, the 
neutralino LSP pair-annihilation is genuinely somewhat suppressed for tan/3 ~ 1 leading 
to too much relic abundance (f2 x ^o > 1) for m 1 / 2 and tuq smaller than in the general case. 
Allowing for larger tan/3 relaxes both lower and upper limits on both TO1/2 an d mo Q and 
thus on the sparticle and Higgs masses. Clearly, the sparticle mass spectroscopy of the 
next generation of colliders will teach us a great deal about our theoretical expectations, in 
particular on the question of b-r unification. 
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